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We  calculate  the  effect  of  local  magnetic  moments  on  the  electron-phonon  coupling  in  BaFe2  As2+5 
using  the  density  functional  perturbation  theory.  We  show  that  the  magnetism  enhances  the  total 
electron-phonon  coupling  by  ~  50%,  up  to  A  <  0.35,  still  not  enough  to  explain  the  high  critical 
temperature,  but  strong  enough  to  have  a  non-negligible  effect  on  superconductivity,  for  instance, 
by  frustrating  the  coupling  with  spin  fluctuations  and  inducing  order  parameter  nodes.  The  en¬ 
hancement  comes  mostly  from  a  renormalization  of  the  electron-phonon  matrix  elements.  We  also 
investigate,  in  the  rigid  band  approximation,  the  effect  of  doping,  and  find  that  A  versus  doping  does 
not  mirror  the  behavior  of  the  density  of  states;  while  the  latter  decreases  upon  electron  doping,  the 
former  does  not,  and  even  increases  slightly. 


PACS  numbers:  63-20. dk,63-20.Kd,74-20.Pq, 74- 70.Xa 

The  simultaneous  presence  of  high  Tc  superconduc¬ 
tivity  and  magnetism  in  the  phase  diagram  of  the  Fe- 
based  superconductors  (FBSC)  suggests  that  magnetism 
plays  an  important  role  in  determining  the  superconduct¬ 
ing  properties.  Phonons  have  been  excluded  early  on  as 
possible  mediators,  on  the  basis  of  first-principles  non¬ 
magnetic  calculations  for  the  undoped  compound,  but 
the  experimental  situation  is  still  far  from  settled  in  this 
regard.  — ^ 

Density  functional  (DFT)  calculations  correctly  repro¬ 
duce  several  properties  of  Fe  pnictides,  such  as  the  mag¬ 
netic  pattern  in  the  parent  compounds  and  the  geometry 
of  the  Fermi  surface.  The  interplay  of  magnetic  and  elas¬ 
tic  properties  is  however  puzzling^  on  one  hand,  exper¬ 
iments  measure  weak  magnetic  moments  (m  ~  0.3  —  1.0 
Hb)  in  the  SDW  state  and  no  long-range  magnetic  or¬ 
der  in  the  superconducting  samples,  while  spin-polarized 
LSDA  calculations  predict  large  magnetic  moments  at 
all  dopings  (to  ~  2.0  /xg).  On  the  other  hand,  the 
same  spin-polarized  LSDA  calculations  predict  equilib¬ 
rium  structures  and  phonon  densities  of  states  that  are 
much  closer  to  the  experiment  than  those  predicted  by 
non-magnetic  calculations.  7~9  A  possible  way  to  recon¬ 
cile  these  two  apparently  conflicting  results  is  an  itinerant 
picture,  in  which  the  Fe  atoms  nevertheless  have  large  lo¬ 
cal  magnetic  moments  that  order  below  the  Neel  temper¬ 
ature  in  the  undoped  compound,  but  survive  locally  at  all 
dopings:10,11.  But  this  point  of  view  immediately  raises 
the  question  of  whether  the  estimate  of  electron-phonon 
(e  —  ph)  matrix  elements  given  by  early  non-magnetic 
(rather  than  paramagnetic,  i.  e.,  with  the  local  Fe  mo¬ 
ment  entirely  suppressed)  DFT  calculations  is  represen¬ 
tative  of  the  actual  compounds^ 

In  this  paper,  we  calculate  from  first  principles  the 
e—ph  coupling  constant  in  antiferromagnetic  (AFM)  and 
non-magnetic  (NM)  BaFe2As2,  using  the  linear  response 
method^  We  confirm  that  magnetism  strongly  affects 
the  phonon  frequencies,  leading  to  a  renormalization  of 


the  modes  that  involve  Fe-As  vibrations, and  also  find 
a  strong  effect  on  the  e  —  ph  matrix  elements,  leading  to 
a  ~  50%  increase  with  respect  to  the  NM  values.  Using  a 
rigid-band  model,  we  show  that  the  e  —  ph  coupling  con¬ 
stant  A  as  a  function  of  doping  does  not  follow  the  density 
of  states.  Finally,  we  estimate  the  e  —  ph  coupling  in  the 
paramagnetic  (PM)  state,  by  combining  the  nonmagnetic 
band  structure  (eigenenergies  and  eigenfunctions)  with 
the  magnetic  phonon  spectra  and  self-consistent  poten¬ 
tials.  For  the  relevant  values  of  doping,  we  estimate  an 
upper  bound  to  the  e  —  ph  coupling  constant  A  =  0.35, 
i.e.  not  high  enough  to  explain  superconductivity,  but 
not  sufficiently  weak  to  be  neglected.  For  instance,  e—ph 
interaction  may  be  one  of  the  factors  responsible  for  ex¬ 
perimentally  observed  gap  anisotropy. 


Magnetic  Order 

m  (/xB) 

A(0)(ev-1) 

Afjcr 

^In  (X) 

Xaa/N{0) 

NM 

0.0 

1.18 

0.18 

194 

0.15 

AFMc 

2.4 

1.36 

0.33 

179 

0.24 

AFMs 

2.6 

0.68 

0.18 

180 

0.26 

PM  1 

* 

1.18 

0.27 

206 

0.23 

PM2 

* 

1.18 

0.27 

195 

0.23 

PM3 

* 

1.18 

0.31 

170 

0.26 

TABLE  I:  Calculated  properties  of  BaFe2As2.  m  is  the  in¬ 
tegral  over  the  cell  of  the  absolute  value  of  the  magnetiza¬ 
tion,  A(0)  is  the  DOS  in  states/(spin  eV  Fe  atom),  Xa<7 
is  the  electron-phonon  coupling  and  ivin  is  the  logarithmic- 
averaged  phonon  frequency.  In  addition  to  the  three  fully 
self-consistent  calculations  we  report  here  three  model  cal¬ 
culations  (see  text):  PM1  utilizes  the  wave  functions,  one- 
electron  energies  and  phonon  frequencies  from  the  NM  calcu¬ 
lations,  and  the  deformation  potentials  from  the  AFMc  cal¬ 
culations.  PM2  and  PM3  use  one-electron  energies  and  wave- 
functions  from  the  NM  calculations,  and  phonon  frequencies 
and  deformation  potentials  from  the  AFMc  and  AFMs  calcu¬ 
lations,  respectively. 
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In  order  to  disentangle  the  structural  and  mag¬ 
netic  effects,  we  used  the  high-temperature  tetragonal 
structure14!15  in  all  our  calculations.  In  order  to  es¬ 
timate  the  effect  of  the  SDW  ordering  vector  on  the 
e  —  ph  properties,  along  with  the  NM  calculations,  we 
have  considered  two  different  AFM  patterns:  the  checker¬ 
board  one  (AFMc),  in  which  the  nearest  neighbor  Fe 
spins  are  antiparallel,  and  the  experimentally  observed 
one  (AFM  stripe,  AFMs),  in  which  spins  are  aligned 
(anti)ferromagnetically  along  the  (■ y)x  edge  of  the  square 
Fe  planes.  The  stabilization  energies  with  respect  to  the 
NM  solutions  are  90  meV  and  130  meV /Fe  atom,  respec¬ 
tively. 

All  calculations  were  performed  in  the  generalized  gra¬ 
dient  approximation^  using  plane-waves  17  and  ultra- 
soft  pseudopotentials^.  We  employed  a  cut-off  of  40 
(480)  Ryd  for  the  wave- functions  (charge  densities) .  The 
electronic  integration  was  performed  over  an  83  k-mesh 
with  a  0.01  Ry  Hermitian-Gaussian  smearing  in  the  NM 
and  AFMc  case,  while  for  the  AFMs  case  we  used  a 
8x4x8  k-mesh  with  a  0.01  Ryd  Hermitian-Gaussian 
smearing.  Finer  grids  (203  and  20  x  10  x  20  )  were 
used  for  evaluating  the  densities  of  states  (DOS)  and 
the  phonon  lincwidths.  Dynamical  matrices  and  e  —  ph 
linewidths  were  calculated  on  43,  NM  and  AFMc  case, 
and  23,  AFMs  case,  uniform  grids  in  q-space.  Phonon 
frequencies  throughout  the  Brillouin  Zone  were  obtained 
by  Fourier  interpolation.  The  (perturbed)  potentials  and 
charge  densities,  as  well  as  the  phonon  frequencies,  were 
calculated  self-consistently  at  zero  doping  (5  =  0);  the  ef¬ 
fect  of  doping  on  the  e  —  ph  coupling  was  then  estimated 
using  the  rigid-band  approximation. 

The  band  structures  at  S  =  0  (not  shown)  agree  with 
previous  calculations]11’19  The  corresponding  magnetic 
moments  and  densities  of  states  (DOS)  at  the  Fermi  level 
are  reported  in  table  The  phonon  dispersion  for  NM, 
AFMc,  and  AFMs  cases  are  shown  in  the  left  panels  of 
Fig.  m  the  remaining  panels  show  the  (partial)  densi¬ 
ties  of  states  and  Eliashberg  functions.  Our  results  agree 
with  previous  calculations  in  the  same  crystallographic 
structures  (adding  orthorhombicity  additionally  changes 
phonon  dispersions,  see  Refs.  @1)).  Magnetism  has  the 
biggest  effect  on  phonon  modes  that  involve  Fe-As  vibra¬ 
tions;  there  is  a  pronounced  softening  of  a  branch  origi¬ 
nally  located  at  ~  25  meV  along  the  T  —  Z  line  in  the  NM 
calculation,  down  to  ~  20  meV  in  both  AFMc  and  AFMs 
calculations.  This  branch  corresponds  to  out-of-plane  vi¬ 
brations  of  the  As  atoms.  The  in-plane  Fe-As  modes 
at  low  and  high  energy  are  sensitive  not  only  to  the  size, 
but  also  to  the  pattern,  of  the  magnetic  moment.  Indeed, 
when  the  order  is  AFMs,  Fe  vibrations  along  the  AFM 
direction  harden,  whereas  those  along  the  FM  direction 
soften,  while  for  the  As  vibrations  it  is  the  opposite. 

The  shift  of  phonon  frequencies  in  the  AFM  calcula¬ 
tions  has  often  been  considered  an  indication  of  an  en¬ 
hanced  e  —  ph  coupling  in  the  magnetic  phase;— but 
no  explicit  calculations  of  the  e  —  ph  coupling  constant 


in  the  magnetic  case  have  been  reported  so  far.  In  this 
work,  we  have  calculated  from  first  principles  the  sarne- 


FIG.  1:  (color  online)  Phonon  properties  of  undoped 
BaFe2As2,  for  different  magnetic  patterns:  non- magnetic 
(NM),  AFM  checkerboard  (AFMc)  and  AFM  stripe  (AFMs). 
From  left  to  right:  the  phonon  dispersions;  the  total  phonon 
densities  of  states,  F(co)  (red,  dashed),  and  Eliashberg  func¬ 
tion,  a2F(oo)  (black,  solid);  partial  contributions  of  Fe  (red, 
dark)  and  As  (green,  light)  ions  to  the  total  F(uj),  projected 
onto  Cartesian  axes  (the  Ba  contribution  is  limited  to  u>  <  10 
meV,  and  not  affected  by  magnetic  order).  The  dispersions 
are  shown  in  the  same  Brillouin  zone  for  the  three  patterns. 
The  k-points  are  selected  so  that  they  are  physically  the  same 
in  all  three  structures;  the  magnetic  structure  makes  the  two 
T  —  M  directions  inequivalent  in  the  AFMs;  the  spins  are 
aligned  AFM(FM)  in  the  x(y)  direction.  In  this  coordinate 
system,  the  points  are:  V  =  (0,0,0);  Z  =  (0,0, 7r/c)  or 
(n/a,  ir/a,  0);  M  =  (n/a,  0, 0);  M'  =  (0, 7r/2a,  0),  where  a 
is  the  length  of  the  shortest  Fe-Fe  bond  and  c  is  the  Fe-Fe 
interlayer  distance. 


3 


spin  component  of  the  Eliashberg  spectral  function: 


(x  F(j(T(uj') 


V,(T 

9  k  n,k+qm 


Na(0)Nk  IVk+ql"  x 

«(e£)<J(^+q)<5(w-0,  (1) 

(k'7n|$T4/£eqi,|k  +  q  a  m)  j  ^[2to^v. 


Here,  Nk  is  the  number  of  k-points  used  in  the  summa¬ 
tion,  N„  (0)  is  the  density  of  states  per  spin  at  the  Fermi 
level,  and  u> q  are  the  phonon  frequencies.  The  e  —  ph.  ma¬ 
trix  element  k+qm  is  defined  by  the  variation  of  the 
self-consistent  crystal  potential  Va  for  the  spin  a  with 
respect  to  a  frozen  phonon  displacement  according  to 
the  phonon  eigenvector  eq„  =  J^Aa  ma  ^/^^e^auqAa- 
Here  uqAa  is  the  Fourier  transform  of  the  a  component  of 
the  phonon  displacement  of  the  atom  A  in  the  unit  cell, 
Ma  is  the  mass  of  atom  A  and  e^a  are  Aa  components 
of  qz'  phonon  eigenvector  normalized  in  the  unit  cell. 
The  first  inverse  moment  of  a2F(to)  gives  the  frequency- 
dependent  e—ph  coupling  constant: 


f-OJ 

\aa(oj)  =  2  dna2Faa(n)/n. 

J  0 


(2) 


The  total  e—ph  coupling  constant  \ac,  =  Xaa (to  =  oo) 
is  0.18,0.33  and  0.18  for  NM,  AFMc  and  AFMs,  respec¬ 
tively. 

The  three  calculations  in  Fig.  |T]have  different  phonon 
spectra  coq,  different  self-consistent  crystal  potentials 
Vfj ,  and  different  one-electron  wave- functions  |kn)  and 
eigenenergies  which  all  determine  the  value  of  A  in 
Eqs.  OS  In  order  to  understand  which  of  these  factors 
dominates  the  e  —  ph  coupling,  we  start  for  simplicity 
from  the  so-called  Hopfield  formula  for  the  e  —  ph  cou¬ 
pling  constant: 


Na(0)D2 
Muj 2 


(3) 


where  D  is  the  deformation  potential,  which  is  a  measure 
of  the  average  e—ph  interaction,  and  tfu2  is  a  character¬ 
istic  force  constant,  which  we  assume  depends  weakly  on 
the  magnetic  order.  We  can  then  use  the  ratio  \aa/Na(fS) 
to  obtain  an  estimate  for  the  average  change  in  the  e—ph 
interaction  due  to  magnetism;  going  from  NM  to  AFMc 
(AFMs),  this  ratio,  and  hence  D2,  increases  by  ~  50%, 
from  0.15  to  0.24  (0.26)  -  see  tabic  [T]  A  comparison  of 
Eq.  [T]  and  Eq.  [3]  shows  that  the  increase  in  D  can  be 
caused  either  by  an  increase  in  the  derivative  of  the  DFT 
potential  Va  with  respect  to  displacement,  or  by  a  differ¬ 
ence  in  the  one-electron  wavefunctions,  used  to  evaluate 
its  average  over  the  Fermi  surface.  In  order  to  disentan¬ 
gle  these  two  effects,  we  performed  a  mixed  calculation 
(PM1),  in  which  we  combined  the  NM  eigenvalues  and 
wavefunctions,  as  well  as  the  phonon  frequencies,  with 
the  AFMc  potential  variations  in  the  Eq.  [T|  In  the  top 
panel  of  Fig.  [U  we  compare  the  resulting  Eliashberg  func¬ 
tion  a2F(ui )  with  the  NM  one;  the  total  phonon  DOS 
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FIG.  2:  ( color  online  :)  The  Eliashberg  functions  a2F(uo) 
for  the  non-magnetic  (NM,  top,  black  solid  line),  anti¬ 
ferromagnetic  (AFM,  middle  and  bottom,  black  solid  lines) 
and  model  paramagnetic  calculations,  as  described  in  the  text, 
(PM1-3,  colored  dashed  lines).  The  corresponding  phonon 
densities  of  states,  in  arbitrary  units,  are  shown  in  each  panel 
by  the  dotted  lines. 


F(to)  is  also  shown  as  dotted  lines.  The  effect  of  us¬ 
ing  the  AFM  potential  is  an  increased  coupling  of  elec¬ 
trons  to  phonons  with  frequencies  oj  ~  20  —  25  meV, 
while  other  modes  are  largely  unaffected.  This  increases 
A  from  0.18  (NM)  to  0.27  (PM1).  The  middle  and  lower 
panel  of  Fig.  [U  show  two  calculations,  PM2  and  PM3, 
which  use  NM  one-electron  wave  functions,  and  AFMc 
and  AFMs  phonon  frequencies  and  crystal  potentials,  re¬ 
spectively.  A  comparison  of  the  a2F(to)  with  the  phonon 
densities  of  states  shows  that,  also  in  these  cases,  there  is 
an  increased  coupling  to  phonons  around  20  meV,  both  in 
AFM  and  PM.  The  A’s,  reported  in  tabic  |H  are  0.27  and 
0.31,  for  PM2  and  PM3  respectively.  From  the  compari¬ 
son  of  PM1  and  PM2,  we  conclude  that  the  effect  of  the 
phonon  frequency  is  negligible,  while  the  ~  10%  spread 
of  values  between  PM2  and  PM3  gives  an  indication  of 
the  effect  of  the  short  range  AFM  correlations  on  A.  Fi¬ 
nally,  we  find  that  the  values  of  \arj /Na(0)  in  the  PM 
calculations  are  in  line  with  the  AFM  ones,  indicating 
that,  at  <5  =  0,  the  states  at  Ep  have  a  comparable,  weak 
coupling  to  phonons,  for  both  NM  and  AFM  orders. 

These  model  calculations,  which  combine  the  AFM 
potentials  and  phonons  and  the  NM  wavefunctions  and 
Fermi  surfaces,  represent  the  best  approximation,  at 
LDA  level,  of  the  real  paramagnetic  (PM)  state  of  su¬ 
perconducting,  doped,  BaFe2As2  which  is  characterized 
by  local,  disordered  magnetic  moments^  In  the  follow¬ 
ing,  we  will  use  them  to  estimate  an  upper  bound  for  the 
e  —  ph  coupling  in  the  doped  BaFe2As2- 
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FIG.  3:  (color  online)  Effect  of  doping  on  the  e  —  ph  coupling 
of  BaFe2As2,  in  the  rigid-band  model.  ±<5  is  the  number  of 
excess  electrons/holes  per  Fe  atom.  The  left  panel  shows  the 
density  of  states  at  Ef,  in  states/spin  eV  Fe  atom.  The  right 
panel  shows  the  same-spin  coupling  constant  A  =  XalT.  Na(0) 
is  the  same  for  NM  and  all  three  (see  text)  PM  calculations, 
as  described  in  text. 


In  Fig.  0  we  show  the  density  of  states  af  the  Fermi 
level  Na{ 0)  {left)  and  Xaa  {right)  as  a  function  of  doping  S 
for  NM,  AFMc,  AFMs,  and  PM  orders,  from  a  rigid-band 
calculation;  5  is  defined  as  the  number  of  excess  electrons 
(holes)  per  Fe  atom.  NM  and  PM  calculations  have  the 
same  A^O),  which  decreases  monotonically  as  a  function 
of  doping,  while  Xaa  is  roughly  symmetric  around  5  =  0. 


This  indicates  that,  in  Eqs.  m  the  effect  of  the  matrix 
element  dominates  over  that  of  the  DOS.  In  the  AFM  cal¬ 
culations,  the  shape  of  the  Xa<7  curve  follows  more  closely 
that  of  Na{ 0),  with  a  minimum  at  5  ~  0.1  —  0.2,  and  a 
maximum  XarJ  =  0.8  for  AFMs  at  5  =  —0.5,  correspond¬ 
ing  to  KFe2As2-  However,  the  use  of  the  rigid-band  ap¬ 
proximation  is  questionable  at  these  high  dopings,  and  in 
the  following  we  limit  our  analysis  to  a  smaller  range  of 
5  =  ±0.25,  where  the  error  in  A  connected  to  the  rigid- 
band  approximation  is  <  20%.— 

We  can  summarize  the  results  of  Fig.  [3j  by  saying 
that,  for  |<5|  <  0.25,  the  upper  bound  for  the  e  —  ph 
coupling  in  paramagnetic,  superconducting  BaFe2As2  + 
25  is  A (j (j  <  0.35.  The  enhancement  in  A  results  from  a 
~  50  %  increase  in  the  e  —  ph  matrix  elements,  due  to 
static  magnetism,  which  is  independent  on  doping,  and 
a  symmetric  increase  of  the  matrix  element  for  5  ^  0, 
which  does  not  follow  the  shape  of  the  electronic  DOS. 

The  value  Xaa  =  0.35,  which  includes  magnetism  and 
doping,  is  almost  a  factor  of  two  larger  than  that  esti¬ 
mated  in  early  non-magnetic  calculations  for  the  undoped 
compounds  A—,  and  close  to  to  a  recent  experimental  es¬ 
timate  of  A  from  a  kink  in  photoemission  spectral 

What  are  the  consequences  of  this  result?  Of  course, 
it  cannot  explain  a  Tc  of  38  K,  which  confirms  that  su¬ 
perconductivity  in  this  and  other  Fe-based  superconduc¬ 
tors  is  most  likely  due  to  electronic  (magnetic)  degrees 
of  freedonii^^  However,  in  Fe  pnictides,  studies  of  the 
superconducting  gap  in  realistic  models  with  AFM  spin 
fluctuations  show  that  solutions  with  and  without  gap 
nodes  are  almost  degenerate,  so  that  even  a  relatively 
low  e—ph  coupling  constant  can  help  select  either  of 
them,  by  enhancing/suppressing  the  pairing  in  the  rele¬ 
vant  channel.24:25 
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